Point defect dynamics in two-dimensional colloidal crystals 
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We study the topological configurations and dynamics of individual point defect vacancies and 
interstitials in a two-dimensional crystal of colloids interacting via a repulsive Yukawa potential. 
Our Brownian dynamics simulations show that the diffusion mechanism for vacancy defects occurs 
in two phases. The defect can glide along the crystal lattice directions, and it can rotate during 
an excited topological transition configuration to assume a different direction for the next period 
of gliding. The results for the vacancy defects are in good agreement with recent experiments. For 
interstitial point defects, which were not studied in the experiments, we find several of the same 
modes of motion as in the vacancy defect case along with two additional diffusion pathways. The 
interstitial defects are more mobile than the vacancy defects due to the more two-dimensional nature 
of the diffusion of the interstitial defects. 
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I. INTRODUCTION 



Topological defects in two-dimensional crystals are rel- 
evant to a number of condensed matter systems includ- 
ing vortices in type-II superconductors Wigner crys- 
tals 0, magnetic bubble arrays [U, atoms on surfaces, 
and dusty plasmas 0]. The creation of topological de- 
fects such as dislocations and disclinations leads to a two 
step melting process and an intermediate hexatic phase 
@)H,0|- Topological defects play a role in the mechanical 
response of the system, such as when a shear is applied to 
the crystal, and also determine how effectively the crystal 
can be pinned to a disordered substrate. More recently, 
there has been growing interest in studying topological 
defects in two-dimensional systems on curved surfaces, 
such as colloidal particles on the surface of a liquid drop 
[1, Q . In addition to providing insight into how defects 
can affect equilibrium and nonequilibrium properties of 
these systems, understanding how individual topological 
defects move would also be valuable for technological ap- 
plications, such as the nanoscale or mesoscale engineering 
of new two-dimensional materials. 

In an effort to understand the dynamics of individual 
point defects in two dimensions, recent experiments were 
conducted in which point defects were artificially created 
in a two-dimensional colloidal suspension by manipulat- 
ing the colloids with optical tweezers [l(| EH G3 • Here, 
a defect is created by removing one or two colloids from 
a perfect triangular lattice, resulting in a mono- or diva- 
cancy. The defect configurations were shown to have a 
lower symmetry than the triangular lattice, and specific 
topological arrangements were characterized by the ar- 
rangement of fivefold and sevenfold coordinated colloids 
around the cores of the mono- and divacancies. The ex- 
periments showed that the trajectories of both types of 
defects followed the major axis directions of the triangu- 
lar colloidal crystal. As a result, at short times the defect 
diffusion had one-dimensional characteristics. 



In this work, we examine the topological configurations 
and dynamics of single point vacancies and interstitial 
colloids in a two-dimensional triangular colloidal crys- 
tal. Only vacancies were considered in the experimental 
work [la El 01, but here 

we compare the behavior of 
vacancies to that of interstitial defects. We find that the 
point defects have numerous stable topological configu- 
rations as well as several frequently appearing excited 
configurations. The diffusive thermal motion of the de- 
fects is aligned with the major axes of the crystal, as also 
seen in the experiments. As the defect moves, it switches 
between gliding configurations, which contain two five- 
fold coordinated particles, and transient configurations, 
which contain three or more fivefold coordinated parti- 
cles. The defect is able to change its gliding orientation 
in the transient configurations. The mobility of intersti- 
tial defects is greater than that of vacancy defects, and 
we show that this is a result of the fact that interstitial 
defects are more likely to undergo reorientations in their 
gliding direction than vacancy defects, giving the diffu- 
sion of the interstitial defects a more two-dimensional 
character. We discuss the implications of our work for 
other systems, such as the effect of the motion of intersti- 
tial and vacancy defects on transport in two-dimensional 
vortex lattices. 



II. SIMULATION METHODS 

We simulate a two dimensional colloidal crystal com- 
posed of N = 1116 ± 1 particles using Brownian dynam- 
ics. The system size is L x — 31ao and L y = 18\/3aoi 
where distances are measured in units of the colloid lat- 
tice constant ao, and where we employ periodic boundary 
conditions in both the x and the y directions. The over- 
damped equation of motion for an individual colloid i 
is 

r^f^fec + ff (1) 
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where the damping coefficient r\ = 1 in simulation 
units. Here the colloid-colloid interaction force is f cc = 
-Y,j^i^iU{rij)rij, where r»j = |r, - is the dis- 
tance between particles located at r; and r^, and fy = 
(r-j — Tj)/rij. We represent the colloid interaction via 
a Coulomb potential that is screened by the presence 
of ions in the liquid phase, giving us the Yukawa form 
Ufrij) = Eoq 2 exp(— Kry )/r,j. The unit of energy is 
Eq = Z* 2 / (Aireeoao) where Z* is the unit of charge, e is 
the solvent dielectric constant, q — 1 is the dimensionlcss 
charge on each colloid, and 1/n = 0.16a is the screening 
length. Our system is in the strongly charged low volume 
fraction limit, allowing us to neglect hydrodynamic inter- 
actions. The thermal force if is introduced as Langevin 
random kicks that obey the equations (if(t)) = and 
(ff(i)fj(t')) = 2r}k B T6 ij 6(t - t'). Temperatures are re- 
ported with respect to the melting temperature /™ for 
the two-dimensional crystal. Throughout this work, we 
consider temperatures in the range 0.37 < fr/f™ ^ 0.62, 
which is high enough to produce observable defect mo- 
tion, but low enough to fall below the temperature at 
which a proliferation of fivefold and sevenfold defect ex- 
citations begins to occur. After initializing the system 
in a triangular lattice, we either add or subtract a sin- 
gle colloid from a location near the center of the sample. 
We then measure the time evolution of the system over 
a long time period of 2 x 10 8 simulation time steps. 

The periodic boundary conditions in our sample pre- 
vent the individual defects from annihilating. Since the 
interaction between topological defects is well known to 
be long ranged, finite size effects are always a concern in a 
periodic system. This is particularly important in studies 
of multiple interacting defects. In our case, we are work- 
ing with only a single defect, which experiences a per- 
fectly symmetrical interaction with its images across the 
periodic boundaries. This tends to minimize the impact 
of the periodic boundaries. To check for finite size effects, 
we tested larger systems and did not observe any changes 
in the diffusive behavior. Additionally, our results for the 
vacancy defects are in good agreement with experiments, 
offering further evidence that we have avoided finite size 
effects. 



III. POINT DEFECT CONFIGURATIONS 

To characterize the defect configuration and identify 
the position of the defect, we use a Voronoi cell con- 
struction performed after the system has relaxed into a 
stationary state from its initial configuration. We can 
identify the coordination number Zi of each colloid by 
counting the sides of the polygons in the Voronoi con- 
struction. In an ideal triangular lattice, all particles are 
sixfold coordinated with Zi = 6. Our system has either a 
missing particle or an extra particle, and the defect core is 
surrounded by a number of z, =/= 6 particles. We identify 
distinct topological configurations of the defect based on 
the Voronoi cell picture, adopting the same notation as 




(a) (b) 




(c) (d) 

FIG. 1: Voronoi cell construction for commonly observed va- 
cancy defect configurations. Colloid positions are indicated 
by dots. The Voronoi cells are colored according to the coor- 
dination number Zi of the colloids: white, Zi — 6; light grey, 
Zi — 7; dark grey, Zi = 5; very dark grey, Zi — 8. Only a 
6ao x 6ao portion of the full system is shown, (a) A twofold 
crushed vacancy (V2a)- (b) A split vacancy (SV) centered 
on a Zi — 8 colloid, (c) A threefold symmetric vacancy (V3). 
(d) A fourfold symmetric excited configuration V4 containing 
eight colloids with Zi 7^ 6. 

used in Refs. [HI HU . All of the possible configurations 
that exist in this system conserve the average sixfold co- 
ordination of the particles, 

N 

J2 z i= 6N - ( 2 ) 

i=l 

In Fig. [Tj we illustrate the most prevalent configurations 
for the vacancy defect. Fig. [Ha) shows that the twofold 
crushed configuration (V2 a ) consists of two nearly paral- 
lel pairs of Zi = 5 and z.; = 7 particles. The split con- 
figuration (SV) for the vacancy, illustrated in Fig. [TJb), 
contains three particles with Zi ^ 6, forming an almost 
straight line with two = 5 particles on opposite sides of 
a zi = 8 particle. The third equilibrium configuration for 
the vacancy, the threefold symmetric configuration (V3) 
shown in Fig.[Tjc), consists of three pairs of particles with 
zi = 5 and Zi = 7 arranged around the outside of a trian- 
gle that is centered on the vacancy. These configurations 
are the same as those observed for the monovacancy in 
the colloidal experiments [TJ EH . Figure DJd) illustrates 
an important excited configuration that we find for the 
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FIG. 2: Voronoi cell construction for commonly observed in- 
terstitial defect configurations. Colloid positions are indicated 
by dots. The Voronoi cells are colored according to the coor- 
dination number Z4 of the colloids: white, z% = 6; light grey, 
Zi — 7; dark grey, Zi = 5. Only a 6ao x 6ao portion of the 
full system is shown, (a) A threefold symmetric interstitial, 
/3. (b) A twofold symmetric interstitial, fa- (c) A disjoint 
twofold symmetric interstitial, fad- (d) A fourfold symmetric 
excited configuration I4. 



vacancy. This fourfold symmetric configuration, termed 
VI, contains eight colloids with Zi 7^ 6. It is short lived 
and plays a role in the mobility of the vacancy, as will be 
described below. 

Figure [5] shows the main configurations for the inter- 
stitial defect. Fig. [DJa) illustrates a threefold symmetric 
interstitial, 73 , composed of a triangular arrangement of 
Zi = 5 and zi = 7 colloids centered around the inter- 
stitial. A twofold symmetric interstitial configuration, 
fa, appears in Fig. HJb). This configuration can split 
to become a disjoint twofold symmetric interstitial, fad, 
shown in Fig. [2jc) . The interstitial defect forms an ex- 
cited configuration that is similar in form to that seen for 
the vacancy defect. Termed a fourfold symmetric excited 
configuration, this interstitial configuration is shown 
in Fig.[2jd). As in the vacancy case, this excited config- 
uration persists only for short times, but as we describe 
in Section V.B, it plays an important role in the motion 
of the defect. 



FIG. 3: Diffusion of a vacancy defect at fr/fr = 0.4. Only 
a 14ao x 14ao portion of the full system is shown. Dark lines: 
trajectory followed by the vacancy over a period of 4.3 x 10 7 
time steps. Light lines: long-time average location of the 
crystalline lattice through which the vacancy is moving. Over 
short times, the vacancy diffuses in the directions of the major 
axes of the crystal. 



IV. DEFECT MOTION AND LOCAL BURGERS 
VECTOR DIRECTIONS 

To study the motion of the two types of defects, we 
must identify the location of each defect. We define the 
defect position to be at the center of the N Zi ^ 6 colloids 
with Zi ^ 6 that are present in the system, 

1 N 

rd =M £r*(l-<^-6)). (3) 

N *<** tt 

With this measure, we can follow the trajectories of a 
defect under thermal diffusion. 



A. Isotropic temperature 

In Fig. [3] we illustrate the trajectory of a vacancy de- 
fect over a period of t = 4.3 x 10 7 simulation time steps 
at a temperature of fr/fr — 0-4. We also indicate the 
time-averaged location of the background triangular lat- 
tice to show that the defect diffusion follows the main 
crystalline lattice directions, [10], [01] and [Tl]. The de- 
fect diffuses by gliding along the lattice directions in a 
one-dimensional random walk, punctuated by occasional 
reorientation transitions. This suggests that at short 
times the diffusive motion has one-dimensional character- 
istics, while for long times the diffusion is isotropic. The 
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FIG. 4: (a) Time series of the orientations ai (upper dark 
line) and 02 (lower dashed line) determined from the two 
local Burgers vectors h[ and b 2 in the gliding vacancy defect 
configurations SV and Vi a for a system with fr/fr = 0.4. 
(b) From the same system, a polar plot of the probability 
P(a) of observing a vacancy defect oriented at an angle a 
with respect to the zero x axis. 



trajectories we observe are very similar to those found 
in the vacancy experiments of Ref. [To| where the defect 
moved along the symmetry directions of the crystal. We 
find that interstitial defects move in a similar manner. 
The reorientation transitions that occur during diffusion 
are mediated by the formation of either a threefold sym- 
metric configuration or one of the excited configurations 
illustrated in Figs. [1] and [2j We describe the reorienta- 
tions in more detail below. 

In order to characterize the motion along the different 
crystalline axis, we require a definition of the orientation 
of the defect as a function of time. The orientation of a 
dislocation is readily obtained using the Burgers vector 
b; however, a point defect has a net Burgers vector of 
b = 0. We can construct local Burgers vectors b^ by 
associating each of the N Zi= $ colloids that have Zi — 5 
with the closest Zi > 6 colloid to form a dislocation, such 
that each Zi — 7 colloid is paired with only one Zi — 5 
colloid, and under the constraint that the total length, 

J2j=i~" 1 1» 1 3 is minimized. 

The twofold symmetric configurations, Via, I2, and l2d, 
along with the split vacancy configuration SV, all have 



only two local Burgers vectors which are close to being 
parallel to each other. For these configurations, which we 
refer to as gliding configurations, we define the orienta- 
tion of the defect to be in the direction perpendicular to 
one of the local Burgers vectors. We also identify the an- 
gle Q.j which each Burgers vector makes with the zero y 
axis. In Fig.[5Ja) we plot the time series of a,\ and a<i for 
a system with a vacancy defect diffusing at /r//™ = 0.4. 
The angles are only defined during the time periods when 
the vacancy is in the configuration Via or SV, but the 
vacancy spends most of the time in one of these two con- 
figurations, which have \a\ — ct2 1 — 7r. Figure H^a) shows 
that «i and «2 remain fixed at a particular angle for ex- 
tended periods of time, punctuated by relatively rapid 
changes to a new angle. A histogram of the a values 
with both a± and «2 combined, plotted in polar coordi- 
nates in Fig. 0|b), indicates that a is correlated with the 
six directions of the crystal lattice. Within our sampling 
error, all six of the angles appear with equal probability. 



B. Anisotropic temperature 

Many two-dimensional systems contain some type of 
anisotropy which could originate in the particle-particle 
interactions or from an underlying weak periodic sub- 
strate modulation [3, [TH, [l(| . Since we found that the 
defects move along the symmetry directions of the crys- 
talline lattice, it may be expected that if some form of 
anisotropy is added, the motion may be more prominent 
along certain directions. A similar system has recently 
been realized experimentally in a colloidal system com- 
posed of superparamagnetic particles [ItJ ■ With the ap- 
plication of an external magnetic field, the interactions 
between the colloids could be made anisotropic. As a re- 
sult, dislocations formed in the system with a preferred 
orientation, resulting in melting along the direction of 
the applied magnetic field. 

We consider the effect of adding an anisotropic tem- 
perature to our system by setting f T / It < 1) where 
f T is fixed at j T j Pt = 0.62. Under these conditions, 
the symmetric diffusion disappears and the defects move 
preferentially in one direction. In Fig. [5] we plot P(a) 
for several different temperature anisotropy ratios. As 
It I St 1S decreased from 1, the value ofP(|a = 7r/2|) in- 
creases dramatically while P{a) for the other four crys- 
tal directions decreases. This indicates that the defect 
spends most of the time with a oriented along the y di- 
rection, meaning that the local Burgers vector is oriented 
along the x direction. Reorientation transitions for the 
defect occur much less frequently than in the isotropic 
case, causing the defect to undergo one-dimensional dif- 
fusion for much longer periods of time, and resulting in a 
highly anisotropic diffusion over time. Such anisotropic 
diffusion is consistent with the results obtained in the 
experiment of Ref. [ItJ ■ 
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FIG. 5: -P(a), the histogram of a values observed during a 
simulation run, for different temperature anisotropics in a sys- 
tem with a vacancy defect. In all cases j^/fr — 0.62. Thick 
continuous line: fr/fr = 0.1, dotted line: fx/fr = 0.3, 
dashed line: fx/fr = 0-5, an d thin continuous line: f T / ' f T = 
0.7. 



V. PROPERTIES OF THE DEFECT DIFFUSION 

We have shown that the defect trajectory follows the 
main crystalline directions. The mechanism of defect 
motion is primarily a gliding process. In the case of a 
vacancy defect, the vacancy can glide in the configura- 
tions SV or V2a , which have parallel local Burgers vectors 
and a well-defined glide direction. The one-dimensional 
gliding motion is interspersed with direction switching 
transitions that occur by means of the V3 configuration 
or the excited Vl configuration. Both of these configu- 
rations have higher symmetry than the gliding configu- 
rations but also have nonparallel local Burgers vectors, 
and hence no well-defined glide direction. Thus, when 
the vacancy enters a transition configuration, the sym- 
metry breaking that defined a glide direction is lost. The 
vacancy spends a relatively short time in the transition 
configuration before reentering one of the glide config- 
urations, which may have the same or a different glide 
orientation than the glide configuration which preceded 
the transition configuration. The high symmetry of the 
transition configurations permit the defect to reorient its 
direction of motion along a new lattice direction. Inter- 
stitial defects undergo similar sequences of gliding and 
reorientation, with the I2 and l2d configurations provid- 
ing the gliding motion, and the -Z3 and 1^ configurations 
permitting the reorientation. 

The vacancy position remains stationary when enter- 
ing a transition configuration such as V3, but the position 
of the interstitial shifts slightly when switching between 
the I2 or l2d and I3 configurations, providing an addi- 
tional mechanism for interstitial motion that does not 
exist for the vacancy defect. It is also possible for in- 
terstitial motion to occur through a 'breathing' process 




Ax, Ay 



FIG. 6: (a) Histograms of the discrete jumps Aa; (solid line) 
and Ay (dashed line) for a vacancy defect, (b) Histograms of 
the discrete jumps Aa; (solid line) and Ay (dashed line) for 
an interstitial defect. In this figure, the jump distances are 
measured in their natural units. Aa; is given in units of ao, 
and Ay is given in units of y/Sao/2 



caused by transitions between the I3 and I 4 configura- 
tions. 

We highlight the difference in the mobility mechanisms 
for the two types of defects by measuring the individ- 
ual jumps of the defects from one lattice position to the 
next. We define Ax, = {v d ^{t + dt) — r dd {t)) ■ x and 
Aj/i = (2/V3){r d , t (t + dt) - r d ,i(t)) ■ y, where we take 
dt = 200 simulation time steps. Histograms of both 
Ax and Ay are plotted in Fig. [Ha) for a vacancy de- 
fect and in Fig. [SJb) for an interstitial defect. The va- 
cancy defect moves only by a glide mechanism. The pro- 
cess of switching between V2 a and SV translates the de- 
fect center by half a lattice constant along one of the 
main crystal axis directions. This produces four non- 
zero values of Ax, ±cos(0)/2 and ±cos(7r/3), and two 
non-zero values for Ay, ± sin(7r/3)/\/3. bmce we are 
plotting the histogram of both Ax and Ay in their nat- 
ural units, which are ao and V3ao/2 respectively, we ob- 
serve peaks at Ax = ±0.25 and ±0.5 and at Ay = ±0.5 
in Fig. [DJa). The same peaks appear for the interstitial 
defect in Fig. Hfb) since the gliding mechanism also op- 
erates in this case, but there are now additional peaks in 
Ax and Ay. These peaks are produced by the other two 
processes that can move the interstitial defect: the switch 
between the I2 or l2 d and I3 configurations, and the 
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FIG. 7: Ar 2 versus time showing the diffusion of the va- 
cancy defect at fr/ fr — 0.62. The dashed line is a linear fit, 
consistent with linear two-dimensional diffusion. 




FIG. 9: (a) Probability to observe different defect configura- 
tions for the vacancy defect. Open circles: P(SV) + P(V2a); 
filled squares: P(Vs); open triangles: P(Vl). (b) Probability 
to observe different defect configurations for the interstitial 
defect. Filled squares: P(h); open circles: P(l2) + P(l2d)', 
open triangles: P(Ii). 
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FIG. 8: Diffusion constant D as a function of temperature 
/t//™- Filled circles: vacancy defect. Open triangles: inter- 
stitial defect. 



breathing motion between the -Z3 and 1^ configurations. 
These two processes move the defect center half a lattice 
constant along directions that lie between the main crys- 
talline directions, at tt/6, tt/2, 57r/6, 77r/6, 3tt/2, and 
ll7r/6. The result is peaks at Ax — ±cos(7r/6)/2 ~ 
±0.43 and peaks at Ay = ±sin(7r/6)/\/3 ~ ±0.29 and 
Ay = ±2 sin(7r/6) /V3 ~ ±0.58. 



A. Differences in defect diffusion for vacancies and 
interstitials 

For short times, the defect trajectory follows the main 
crystalline directions, as was shown in Fig. [3] On the 
longer time scales, however, we find a linear diffusive be- 
havior when we measure the distance that the defect has 



traveled from its original location as a function of time, 
Ar 2 = \r d (t) - r d (0)| 2 . This quantity is plotted in Fig. [3 
for a vacancy defect at faff™ = 0.4. We find that Ar 2 
increases linearly with time, as indicated by the fit in 
the figure. Linear diffusion of monovacancies was also 
observed in the experiments of Ref. . 

To compare the mobility of interstitial and vacancy 
defect, we measure the diffusion constant D, given by 

D _ ^ \rd(t + dt)-r d (t)\ ^ (4) 

with dt — 200 simulation time steps, at a series of tem- 
peratures for each defect type. The results are plotted 
in Fig. [8l The diffusion constant increases with temper- 
ature in each case. Importantly, the interstitial deject 
diffuses significantly faster than the vacancy defect at all 
temperatures. 

B. Probability of different defect configurations 

To help account for the faster diffusion of the intersti- 
tial defect compared to the vacancy defect, we compare 
the time the defects spend in the different topological 
configurations. In Fig. Uta) we plot the probability P of 
observing each configuration as a function of tempera- 
ture fr/ fr f° r the vacancy defect, and compare this to 
Fig. [9^b) which shows the probability of observing inter- 
stitial defect configurations versus temperature. 

Figure [9ja) shows that the vacancy defect spends most 
of the time in the split (SV) and the twofold crushed 
(^2a) configurations, which both allow for gliding motion. 
The excited configuration V[ occurs very infrequently 
and does not contribute significantly to the diffusive pro- 
cess. As the temperature increases, the transition states 
V3 and V{ occur with greater probability, indicating that 
the defect can reorient more easily. A greater number of 
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defect reorientations combined with higher thermal exci- 
tation leads to a higher diffusion constant at higher tem- 
perature. Since the gliding motion is one-dimensional, 
the defect can only move away from its starting position 
as rapidly as a one-dimensional random walk if no re- 
orientation occurs. Reorientation of the glide direction 
produces a more two-dimensional diffusion which allows 
the defect to move away more rapidly from its starting 
position. 

The interstitial defect spends much more time in the 
transition configuration 1% than in any other configura- 
tion, as indicated in Fig. [DJb). The excited I' A configu- 
ration for the interstitial defect also appears with much 
higher probability than the excited V( configuration for 
the vacancy defect case. Overall, this produces a much 
higher probability for reorientation of the direction of 
motion of the interstitial defect, and makes its diffusion 
more two-dimensional than that of the vacancy defect at 
all temperatures considered here. The result is a greater 
mobility of the interstitial with respect to the vacancy. 



VI. DISCUSSION 

An open question is how general our results are for 
defect motion in other two-dimensional systems where 
the particles have different types of interaction poten- 
tials. Although it is beyond the scope of this paper 
to address this issue in detail, we believe that our re- 
sults should hold for any two-dimensional system with 
relatively soft particle-particle interactions in which the 
ground state is a triangular lattice. Very short range 
potentials, such as hard sphere interactions, are likely 
to produce significantly different behavior. It is possi- 
ble for two-dimensional systems to have some other type 
of ordered ground state, such as a square or rectangular 
lattice. These lattices have different available gliding di- 
rections compared to the triangular lattice, so there could 
be either fewer or more modes of motion available to the 
defects. This would be an interesting issue to explore in 
the future. 

Our results could have some implications for the mo- 
bility of defects in other two-dimensional systems where 
the number of vacancies or interstitials can be carefully 
controlled, such as superconducting vortices [HI, EH or 
colloids [2(| interacting with periodic substrates. At par- 
ticle densities where there is one particle per substrate 
minimum, termed a matching density, the system is free 
of topological defects, while vacancies appear for slightly 
lower particle densities and interstitial defects form at 
slightly higher particle densities. If the interstitial defects 
remain more mobile than vacancy defects in the presence 
of a substrate, the interstitials should depin more readily 
than the vacancy defects. In the case of superconducting 
vortices, this would imply that the critical current drops 
off more rapidly above the matching density (when in- 
terstitials are present) than below the matching density 
(when vacancies are present). This is agreement with 



earlier observations [19(. 

Vacancy and interstitial defect motion has also been 
studied in three-dimensional systems in the context of 
atomic crystals. Such systems are not only massive, mak- 
ing them sensitive to vibrational modes, but also can 
form a variety of different lattice structures depending 
on the details of the atomic interactions. A convenient 
method for producing isolated vacancies or interstitials 
in an atomic crystal is through ion damage processes. 
We note that in this context, unusually high mobility of 
self-interstitial atoms (SIAs) has been observed and at- 
tributed to the high likelihood of rotational motions for 
the dumbbell atomic configurations associated with the 
SIA [21( 1 . Although the dimensionality of the resulting 
motion is different, this resembles our observation that 
interstitial defects are more mobile than vacancy defects 
due to the greater likelihood of defect reorientation. 

Another issue is the relative stability of interstitial and 
vacancy defects. Due to the geometry of our system, 
which has periodic boundary conditions, both the inter- 
stitial and vacancy defects are completely stable since 
there is no free surface where the defects can annihilate. 
In a real crystal, the more mobile species of defect will 
reach the edge more rapidly and annihilate. Thus, one 
implication of the faster mobility of interstitial defects is 
that a crystal that forms with an initial population of va- 
cancy and interstitial defects will, over time, experience a 
greater loss of interstitial defects to the crystal edge com- 
pared to vacancies, particularly if there is a strain field 
that helps drive the defects to the edge. Both interstitials 
and vacancies can still annihilate within the crystal, but 
the difference in edge annihilation could result in a dom- 
inant population of the slower moving vacancy defects in 
the relaxed crystal. 



VII. CONCLUSIONS 

In summary, we have studied the topological configura- 
tions and dynamics of individual vacancy and interstitial 
defects in a triangular two dimensional colloidal crystal. 
We use a Voronoi cell construction to characterize the 
different topological configurations assumed by the de- 
fects. For vacancy defects, we find the same configura- 
tions that were observed in recent experiments. We show 
that interstitial defects, which were not studied in the 
experiment, appear in distinct configurations and have 
a significantly higher mobility than the vacancy defects. 
The mobility is affected by short lived excited config- 
urations which have four or more fivefold or sevenfold 
coordinated colloids surrounding the defect. The defect 
diffusion process is one dimensional at short time scales 
when the defects glide along the symmetry directions of 
the crystal. Periodically, the defects enter a transition 
configuration which has higher symmetry than the glide 
configuration. This permits the defect to reorient its glide 
direction when it returns to a glide configuration after a 
short period of time. As a result of these reorientations, 
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linear two-dimensional diffusion occurs in the long time 
limit. Application of an anisotropic temperature causes 
the diffusion process to favor one of the permitted gliding 
directions. We describe the mechanisms of defect trans- 
lation in detail, and show that the vacancy defect has one 
mode of motion, but the interstitial defect has three. It is 
these additional modes of motion that give the interstitial 
defect a has a much higher probability of undergoing a re- 
orientation transition than the vacancy defect at a given 
temperature, causing the interstitial diffusion to be more 
two-dimensional in character and therefore faster. 
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